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SUMMARY 


The description of a transient analysis program for computing structural 
responses to input base accelerations is presented. A "hybrid" modal formula- 
tion is used and a procedure is demonstrated for generating and writing all 
modal input data on user tapes via NASTRAN. Use of several new Level 15 
modules is illustrated along with a problem associated with reading the post- 
processor program input from a user tape. An example application of the pro- 
gram is presented for the analysis of a spacecraft subjected to accelerations 
initiated by thrust transients. Experience with the program has indicated it 
to be very efficient and economical because of its simplicity and small central 
memory storage requirements. 


INTRODUCTION 


Design loads in aerospace subassemblies or components are often specified 
in terms of induced acceleration at the mounting interface. This concept has 
been traditionally used to qualify aerospace hardware by subjecting it to 
prescribed input accelerations on a vibration exciter. Transient analyses of 
subassemblies for prescribed acceleration inputs at the interface are, there- 
fore, valuable for designing and augmenting vibration tests, and for computing 
design loads where vibration tests are not practical. 

The transient analysis in the current level of NASTRAN (Level 15.1) does 
not directly provide for input acceleration forcing functions. By using the 
artifice of placing a large mass (with respect to the total system mass) at the 
desired acceleration input point, an input force equal to the total mass times 
the prescribed acceleration will approximate an acceleration input. Theoreti- 
cally, as the fictitious added mass becomes Infinite, the answer becomes exact. 
But as the mass becomes large, the mass matrix tends to become ill conditioned. 
Experience has indicated the "fictitious mass" approach is not desirable. 

In addition, the current NASTRAN transient analysis allows for initial 
conditions only in the direct formulation. The modal formulation, which is 
generally faster and more economical to run, assumes zero Initial conditions. 

In transient analyses of prestrained structures (such as a missile just prior 
to a burnout transient), the initial conditions become very important in pre- 
dicting the magnitude of the structural loading. 
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The purpose of this paper is to describe a transient analysis program 
which has been developed Lo circumvent the abovementioned NASTRAN limitations. 
This program employs a modal formulation and allows for nonzero initial condi- 
tions. It is assumed that the modal input data have been generated and written 
on user tapes by NASTRAN. Hence, the program is termed a postprocessor program. 


A complete derivation of the program theory is presented along with a 
detailed discussion on the generation and reading of NASTRAN user tapes. This 
exercise demonstrates the versatility of several of the new modules added to 
Level 15.1 as well as some of the limitations of user tapes. Finally, an example 
application of the program to the transient analysis of a spacecraft is presented. 
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SYMBOLS 

rigid body transformation matrix for £, set 

internal members load vector 

identity matrix 

stiffness matrix 

mass matrix 

generalized mass matrix (equation BIO) 

coupled flexible body, rigid body mass matrix (equation 5) 
matrix of modal element force vectors 

expansion of rigid body transformation matrix to g set 
vector of displacement components 

time derivation of modal coordinate vector (equation 10) 
critical viscous modal damping coefficient matrix 

eigenvalue matrix 

matrix of modal eigenvectors 
vector of modal coordinates 


Subscripts: (See Appendix A) 

ex subset of total members in structure 

Y subset of g set 
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Notation: 


[ 1 

[ ] 
N 


square or rectangular matrix 


transpose of matrix 
diagonal matrix 
column matrix 


PROGRAM THEORY 



In this section, the theoretical basis for a program to compute the 
transient response of a structure to acceleration forcing functions is given. 
The equations of motion are developed in terms of a "hybrid" modal formulation 
and reduced to a form which makes maximum use of NASTRAN generated eigenvalue 
data. Numerical solutions to the resulting equations are discussed along with 
treatment of the initial conditions. Finally, equations are presented for con- 
verting the modal response data into transient member loads and grid point 
accelerations . 


In the derivations an attempt has been made to generally utilize the nota- 
tion presented in the NASTRAN Manuals (ref ere ces 1, 2, and 3) for ease of 
reading and implementation of the resulting t lations. In particular, the set 
notation of Appendix A, which is taken from Section 1.7.3 of reference 3, is 
used throughout; although, the r set has a somewhat different meaning herein. 
This difference will become apparent in the course of the derivation. 


Equations of Motion 


Assuming no external loads are acting on the grid points of a structural 
system, the undamped equations of motion for displacement set ("a) become 


HQ * HQ ' 0 


(i) 


where 


C“J and are the reduced mass and stiffness matrices, respec- 

tively (see Section 3.5 of reference 2). It is assumed that the system described 
by equation (1) is not completely constrained against rigid body motions (i.e. , 
it can have from 1 to 6 rigid body degrees of freedom). Equation (1) may be 
partitioned as follows: 


Mfl 0 I Mn 

h i-r-i - -i 

M i ! M rr 



- 0 


( 2 ) 
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if con- 


where by definition, the subset {u } of the displacement vector {u }, 

r 3. 

strained, would be just sufficient to eliminate rigid body motion without intro- 
ducing redundant constraints. Selection of the subset {u f } is arbitrary and 

for the present analysis it is chosen to correspond to the input acceleration 
degrees of freedom (a.d.o.f.), and it is specified on a NASTRAN "SUPORT" Bulk 
Data Card. It should be noted that by using the r set for input accelera- 
tions, the a.d.o.f. are restricted to a maximum of six. This restriction is not 
a major limitation since the base of many components can be assumed to be rigidly 
constrained to a plane. The redundant points in the base can thus be assumed 
to be rigidly attached to a single acceleration input point. 

The mathematical problem at hand is to determine the transient response 
of the {u^ } subset to prescribed (u^,) inputs. A solution using a modal formula- 
tion is presented in the following. This approach allows a significant reduction 
in size of the problem with little loss in accuracy by truncating the number of 
modes included in the solution. 

Modal Coordinate Transformation 

The following "hybrid" transformation between modal coordinates (£) and 
physical coordinates (u) is introduced: 



is the matrix of eigenvectors of the structure with the (u r ) coordinates 

constrained to zero (see Appendix B) ; QQ is the identity matrix; and U i } 
is the vector of flexible body modal coordinates. 

Modal Equations of Motion 

Substituting equation (3) into equation (2) ; premultiplying by the trans- 
pose of the transformation matrix, and using equations (B5) , (B6) , (BIO), and 
(Bll) of Appendix B leads to the following 
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' -•* CM 


and 


[”tr] ■ [4 M uV + M L D lr + D L M lr +M rr] 


( 6 ) 


The upper partition matrix of equation (4) yields 

- 1 . 


{^i> + {C i > - ^M jL J C M irl^ u r ^ 


(7) 


Adding viscous modal damping to equation (7) yields the desired equation 
of motion of the system as 


<«1> + (? 6 l“i3<e i } + C“i 2 J {5 i ) ■ -M' 1 [ M ir] { \ } <»> 


where 8^ is the critical viscous damping ratio for the ith mode. 

With the exception of the 6^ values all of the other coefficient values 

are easily obtained as output quantities from a NASTRAN normal mode analysis 
(Rigid Format 3) . 


Method of Solution 

The method of solution used in the program to solve the equations is a 
standard fourth order Runge-Kutta numerical integration routine with variable 
step size error control. Use of this subroutine required reduction of equation 
(8) to first order and generation of the initial conditions in terms of modal 
coordinates. These procedures are discussed in the following sections. 

Reduction to first order equations . - Integration via the Runge-Kutta Sub- 
routine requires the system equations to be a set of first order differential 
equations of the form 

* * f j( y l» y 2» * * * y n^ j - 1, 2, ... n ( 9 ) 

Equation (8) can be transformed to the form of (9) by introducing the auxiliary 
variable {v^} where 

Ui> - (v i ) (10) 

Using equation (10 )» equation (8) than leads to 
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( 11 ) 


<V ■ - C"l 2 J ( «i ! - C 1/M li][ M ir] ,ii r> 


Equations (10) and (11) are now a set of equations in the form of equation (9) 
(as required) and are integrated simultaneously. 


Initial conditions . - If the initial conditions are known for each of the 
(u } coordinates, then the modal initial conditions can be determined by pre- 

3 . 

multiplying equation (3) by the matrix 



it 


rr 


* 


using equations (B5) , (B6) , and (Bll) ; and solving for {£.} to obtain 



By taking the time derivative of both sides of equation (12), the initial condi- 
tions for {|. } can also be determined in terms of {u.} and {u } values, 
x k r 

Since the initial conditions are not generally known in terms of the {u fl } 

coordinates, a different approach was taken for the present program. This 
program assumes 

{^( 0 )} - 0 ( 13 ) 

and 

{v i (0)} - o (14) 

that is, the structure is assumed to be initially in a steady state deformed 
position. The initial conditions are then computed from equation (ll) as 

V°» • -[l' ii ll“lj[ M ir] Ia r (0)) (15) . 

Equations (14) and (15) thus yield the necessary, Initial conditions 
for numerical integration of equations (10) r.nd (11) . 
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Computed Response Data 


Having obtained the transient response of the modal coordinates, it is 
desirable to transform these variables into transient member loads and accelera- 
tions in terms of the physical {u^} coordinates. 

The member loads are given by superposition as 


lf 0 > - ( 16 ) 

where M is a matrix of modal element force vectors for an arbitrary subset 

a of the total members in the system and {F^} is the total load vector corres- 
ponding to the subset a. The modal member load vectors are obtainable as stan- 
dard output from a NASTRAN Normal Mode Analysis. 

Using equation (10), the grid point accelerations are computed from the 
second time derivative of the upper partition of equation (3) to be 



[StfV + Dvpr 1 


(17) 


where y is an arbitrary subset of the g set and |RB gr J is a merger of 

the rigid body transformation matrix with the o , s , and m sets. 

The matrices and are generated by NASTRAN and the subject 

postprocessor program selects the subset y to be picked up for use in equation 
(17). The quantity {v^} is given by equation (11) and (u^; is a gi ve n input 
vector. 

If the quantities (u } and {u^} are also desired, equation (17) can be 

reduced to a set of first order equations to be integrated simultaneously with 
equations (10) and (11). For the present program, however, this was not done. 


NASTRAN GENERATED INPUT DATA 


The rigxd body and flexible body modal data necessary for solution of the 
foregoing equations are easily generated by a NASTRAN Normal Mode Analysis and 
written out on user tapes. Hence, the present postprocessor program is designed 
to read the majority of its input directly from the NASTRAN user tapes. 

A DMAP alter package is required to generate part of the NASTRAN data and 
to write the user tapes. Also, interrogation of the user tapes to read the 
data requires special considerations. Both of these aspects are discussed 
in the following sections. 
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DMAP Alter Package for Normal Mode Analysis 

A listing of the DMAP alter package is given in Appendix C. A brief dis- 
cussion and explanation of the most significant statement is given in the 
following tabulation. Note the use of the new Level 15 modules VEC, UMERGE, 
and 0UTPUT2. (See sections 3,4, 5.2, and 5.3 of ref. 1 and section 3.5 of ref. 2) 


DMAP Alter 

Statement No. Function 


2 and 3 


4 through 20 


6 


7 and 8 


12 and 13 


17 


18 and 19 


21 


Merge an r x r null matrix with the 
D^ r matrix to create a pseudo a x r 

size rigid body modal matrix 

ES* KJ 

Merge the omitted coordinates, single 
point constraint coordinates, and the 
multipoint constraint coordinates, if 
present, with to obtain £rb^ 

Recovery of omitted coordinates 

■ [ G o.][*V] 


Merging of omitted coordinates 


RB 

ar 

-> 

[rB 1 

RB “ 


fr 

or 


L J 


Merging of SPC constraints assuming they 
Are zero 



Recovery of dependent MPC coordinates 

(?V] - OvTnr] 

Merging of dependent MPC coordinates 



mnt. [RBj 
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DMAP Alter 
Statement No. 


Function 


22 and 23 
24 and 23 


26 and 27 


Write on a user ta P e 

Change the eigenvalue problem from "a" 
size to "X." 9ize in agreement with equa- 
tion ( 59 ) 

Change the checkpointed modal matrix 
from "a" size to "X." size 


29 


30 and 31 
32 




Merging of the 
to zero with the " 


£i 

o“ 



coordinates constrained 
X," size modal matrix 


O J ' [*J T [ M U D tr + 


Prints [M Ir ] 


33 


Writes on a user tape 


35 and 36 Write the modal deflections (0PH1G) , 

the modal SPC forces (0QG1) , and the 
modal element forces (0EF1) on a user 
tape 


Interrogation of User Tapes 

When using NASTRAN user tapes for input to postprocessor programs, the 
analyst must read the tape and selectively extract the required input from the 
totality of data present. This task requires either a prior knowledge of the 
format used in writing the tape or interrogation of the tape to see how it is 
written. 


Unfortunately, the user tapes generated by NASTRAN are written in unfor- 
mated binary (i.e. with a mixture of integer, floating point, and alphanumeric 
formats). In addition, some of the data is packed (i.e., zeros omitted). This 
randomness eliminates a prior knowledge of the format. 

Interrogation of the tapes using standard tape dump routines is also some- 
what futile since these programs read and print all data in a single format. 

To illustrate this problem, 19 records of a typical user tape, written to an E 
format, are listed in figure 1. Obviously, much of the data given is meaning- 
less. 
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This interrogation problem was circumvented for the present program by 
writing a special tape dump program for the Langley Research Center CDC 6000 
series computer to reed and print the mixed format. The program logic was 
patterned after the NASTRAN module TABPRT and a listing is given in Appendix 
D. Using the Appendix D program, the same tape used to generate figure 1 was 
again read and the results are given in figure 2. From this improved interro- 
gation the analyst can easily find where desired data are located and adjust 
the read statements in the postprocessor program accordingly. 

From the foregoing discussion, it is apparent that a postprocessor pro- 
gram must be dynamic. That is, the input read statements must be continually 
changed to fit each new problem after interrogation of the user tape. 

EXAMPLE APPLICATION 

The subject transient analysis program was developed in support of the 
Viking Project, which has a mission to soft-land a scientific payload on the 
surface of Mars in 1976. In particular, this program was intended to provide 
transient loads and accelerations in the Viking Dynamic Simulator (VDS) shown 
in figure 3 for input acceleration transients at the base of the Centaur truss 
adaptor. The VDS is a dummy spacecraft, which is dynamically similar to the 
actual Viking spacecraft, and will be flown on a proof (or test) flight of 
a new Titan D-1T Centaur launch vehicle configuration in 1974. This launch 
vehicle will be used for the Viking mission and is shown in figure 4 along with 
the Viking spacecraft. 

Several discrete transient events Induce high loads into the VDS with the 
more prominent of these being Titan Stage 0 Ignition, Titan Stage 1 Shutdown, 
and the Centaur Main Engine Cutoffs. All of these events were analyzed in 
detail using the subject program and some of the results from the Titan 
Stage 0 Ignition were selected as a typical illustration of input and output 
data. 


Six degree-of- freedom acceleration Inputs into the base of the VDS were 
determined analytically from a transient loads analysis of the actual Viking 
configuration as depicted in figure 5* Input to the transient analysis of the 
actual Viking configuration was based on measured force transients from previous 
Titan launches. The base of the VDS was constrained to a plane and the six com- 
ponents of acceleration were input at a single grid point in the center of the 
base. 


A typical set of input translational components of acceleration are 
shown In figure 6. The longitudinal (or Z) component is seen to be the most 

2 

significant and it starts at lg (9.81 m/sec ) to represent the initial gravity 
load on the vehicle resting on the launch pad. This gravity loading causes an 
initial condition on the modal coordinates {£^(0)} as indicated in equation (16). 

The near sinusoidal oscillation of the Z component after 0.5 second is attri- 
buted to excitation of a longitudinal mode of the vehicle by the initial 
thrust transient. 
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Using che input acceleration (see figure 6) and the modal data (for 18 
modes) from the NASTRAN analysis, selected loads for VDS members and 
accelerations were computed using a viscous damping model which is a function 
of the modal frequencies. A typical load-versus-time response is shown for a 
member of the Viking Spacecraft Adaptor in figure 7. Similarly, the transla- 
tional terms of the acceleration computed for the top of the VDS, are shown 
in figure 8. 

A typical computer run for the example problem on the Langley Research 
Center CDC 6600 Computer required a storage of 53 000g and 200 CPU seconds 
fur execution (Including time for generation of 31 output plots) . A comparable 
NASTRAN transient analysis would require at least 120% increase in storage, and 
considerable increase in CPU time and calls to the operating system. Postpro- 
cessor programs thus are seen to offer potential economic savings in addition 
to special purpose capability. 


CONCLUDING REMARKS 


The transient analysis program described herein yields a simple, convenient, 
and economical approach for treating input accelerations and modal initial condi- 
tions. Other than the limitation of six on the maximum number of input accelera- 
tion components, the program is applicable to a broad spectrum of structural 
applications. 

The fact that such a postprocessor program could be simply written to inter- 
face with NASTRAN demonstrates the expanded utility of NASTRAN via the new Level 
15 utility modules and user tape option. Tailor-made programs such as the pre- 
sent one can be designed to be very efficient in comparison to NASTRAN. Thus, 
the authors would encourage further additions and refinement of postprocessor 
convenience modules rather than expanded capability and complexity of NASTRAN. 

In particular the formats for writing user tapes, so that they may be easily 
read by postprocessor programs, should be given prime consideration. 
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APPENDIX A 


THE NESTED VECTOR SET CONCEPT USED TO REPRESENT COMPONENTS OF DISPLACEMENT 


In constructing the matrices used In the Displacement Approach, each row 
and/or column of a matrix is associated closely with a grid point, a scalar 
point, or an extra point. Every grid point has 6 degrees of freedom associated 
with it, and hence 6 rows and/or columns of the matrix. Scalar and extra points 
only have one degree of freedom. At each point (grid, scalar, extra) these 
degrees of freedom can be further classified into subsets, depending on the 
constraints or handling required for particular degrees of freedom. (For example, 
in a two-dimensional problem all z degrees of freedom are constrained and 
hence belong to the s (single-point constraint) set.) Each degree of freedom 
can be considered as a "point," and the entire model is the collection of these 
one-dimensional points. 

Nearly all of the matrix operations in displacement analysis are concerned 
with partitioning, merging, and transforming matrix arrays from one subset 
of displacement components to another. All the components of displacement of 
a given type (such as all points constrained by single-point constraints) form 
a vector set that is distinguished by a subscript from other sets. A given 
component of displacement can belong to several vector sets. The mutually 
exclusive vector sets, the sum of whose members are the set of all physical 
components of displacements, are as follows: 


points eliminated by multipoint constraints 
points eliminated by single- point constraints 
points omitted by structural matrix partitioning 

points to which determinate reactions are applied in static analysis. 


u. the remaining structural pointa used in static analysis (points left 
1 over) 

u extra degrees of freedom introduced in dynamic analysis to describe 
e control systems 

The vector sets obtained by coobinlng two or more of the above seta are 
(+ sign indicates the union of two sets) 

u ■ u + u» , the set used in real eigenvalue analysis 

A T * 

u . » u + u , the set used in dynamic analysis by the direct method 
dll 


t 
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-5 


u. ■ u + u , uncor.:»r< kilned (f! ree) structural points 
f a o 

u ■ u r + u , a) i structural points not constrained by multipoint constraints 
n r s 

u * u + >.» . .1 ! : structural (arid) points including scalar points 
g n r r 

°P ** U g + physical points 

In dynamic analysis, additional vector sets are obtained by a modal trans- 
formation derived from real eigenvalue analysis of the set u . These are 

s 

rigid body (zero frequency) modal coordinates 

Cf finite frequency modal coordinates 


£. ■ £ + the set of all modal coordinates 
i o > f* 

One vector set is defined that combines physical and modal coordinates. 

The set is + u fi , the set used in dynamic analysis by the modal method. 

The nesting of vector sets is depicted by the following diagram: 



The data block USET (USETD in dynamics) is central to this set classifica- 
tion, Each word of USET corresponds to a degree of freedom in the problem. 

Each set is assigned a bit in the word. If a degree of freedom belongs to a 
given set, the corresponding bit is on. Every degree of freedom can then be 
classified by analysis of USET. The coamon block /BITP0S/ relates the sets 
to bit numbers. 
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APPENDIX B 


MODAL PROPERTIES 


In this section, several identities relating to both the rigid body modes 
and the flexible body modes are presented. Although these identities are 
perhaps familiar, they are included herein for completeness an. continuity 
of notation. 

Rigid Bocl* Modal Properties 

For periodic motion of frequency w, equation (2) reduces to the eigenvalue 
equation 


r i 

v 1 If 

U « K lr 

2 

— ■ 

’ M U ! M ir 


t" * ” ” 

k: <k 
L ** i rrj 

u 

- 0) 

mT 1 M 

L * r i rr 



0 


(Bl) 



The solution of equation (Bl) yields the natural frequencies and the corres- 
ponding natural modes of the system. For the rigid body modes corresponding 
to u - 0, equation (FI) reduces to 


(B2) 


Since the rigid body mode matrix p £r ] relates the rigid body motions 
{ u ^} ^ in terms of ^ u r ^ ri g » the following transformation may be written: 

M Vrl,. <63> 

where [ij is the identity matrix. Using equation (B3), equation (B2) gives 


>w j hi 

r* 

hr_ 

5r , K rr 

1 


'“r'rl,. 


(BA) 
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For arbitrary {u } . displacements, it follows from the partitions of (B4) 
r rig* 

that 

M?*r] + for] ■ 0 CBS) 


and 

Kr JfoJ + ILr] ' 0 (B6) 


Equations (B5) and (B6) thus yield two important identities relating the 
rigid body modes and the partitions of the stiffness matrix. 

It should also be noted that solving equation (B5) for the rigid body 
mode matrix yields 


■ -CvTCv] < b 

which is consistent with equation (41) in Section 3.5 of reference 2 and is 
the equation used in NASTRAN Rigid Format 3 to compute the rigid body mode 
matrix. 


Flexible Body Modal Properties 

By definition of the {u r } degrees-of- freedom, introduction of the 
constraint 

{u r ) - 0 (B8) 

eliminates rigid body motion leaving only flexible body motion. Using equation 
(B8), the upper partition of equation (Bl) yields the following eigenvalue 
equation for the flexible body modes: 

" “ M AbJW* " 0 (B9) 

The modal matrix £$£^3 t ' ie * eigenvectors of equation (B9) is shown 
in reference 4 to satisfy the following orthogonality relationships! 


PAW ■ 

‘ [ViS 


The above equations provide the foundation for modal formulation. 
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APPENDIX C 


DMAP ALTER PACKAGE 


01 

02 

03 

0L 

05 

06 
OT 
08 

09 

10 

11 

12 
13 
lU 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2U 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3fc 

35 

36 

37 


ALTER bl 

vec uset/va/c*n*a/c*n*l/c*n*r % 

MERGE UM*t» * *VA/R6A/C*N*0/C*N»2/C*N»S 4 
EQUI V HBA jRBF/OMIT 4 
CONU LBOMtOMIT S 

MPYAU GO,HBA,/RoO/C t N*0/C*N t l/C t NtO/C*N,l 4 
VEC USET/VF/C *N*F /C*N* ,",/C*N*0 4 

MERGE RBA*RBO****VF/RBF/c*N*0/C*N*2/C*N*2 4 
LABEL LBOM $ 

EQUIV RBFtRBN/SlNGLE 4 

CONU LBGStSINGLE 4 

VEC USET/VN/C*N*N/C*N*f/C*N*S 4 

MERGE RBF ♦ ♦ t ♦ * VN/RBN/C *N* 0/C *N*2/C*N*2 4 

LABEL LBSS 4 

EOUIV RBN*RBG/MPCF2 4 

COND lBMM,MPCF2 4 

MPYAU GM*RBN*/RBM/C*N*0/C*N*1/C*N*0/C*N*1 4 
VEC USET/VG/C*N*G/C*N*K/C*N*M % 

MERGE RBN» RBMt * * * VG/RbG/r. *N*0/C*N*2/C*N»2 4 

LABEL LBMM 4 

MATGPR GPE*USET*4IL*RBG//c*N*R/C.N*G 4 
OUTPU72* *«**//C*N*-l/C*N*ll 4 
0UTPUT2* RBG****//C*N*0 /C*k* 11 4 
ALTER 89*b9 

READ KLLfMLL* ♦ *EE0* *CASFCC/LAMA*PH1L*MI*0EIGS>/ 

ALTER 91*91 

CHKPNT LAMA*Ph1L*M1*OE1GS 4 

ALTER 93 

UMEROE USET*PHlL*/PHIA/V*N*MAJOR=A/V*N*SUBO*L/V*N*SUBl=R s 
MPYAU MLL*OM*MLR/TMP/C*N,0/C*N*1/C»N»1/C*N*1 4 
MPYAU PHIl*TMP*/M1R/C*N*1/C*N*1/C*N*0/C*N*1 $ 

MATPHN MIR****// » 

OUTPUTS* M1R*»»*//C*N*0/C*n*11 4 
ALTER 105 

OUTPUTS* OPHIG*OQG1*OEF1**//C*N*0/ *N*11 4 
OUTPUTS* ♦ * • *//C*N*-9/C*N* \ 1 4 

enualter 



»f»r>o»<von*o on 





APPENDIX D 

KUl.TI-l’OUMAT TAPE DUMP PROGRAM 

PROGRAM A3930 (INPUT. 0UTPuT.TAPE5=INPUT,TAP£6=0UTPUT*TAP£l> 


c 

THIS PROGRAM READS any TAPE WITH VARIA3LE LENGTH BINARY RECORDS C 

SLOWEST TAPE MUST BE LADf L ED AS TaPEI C 

7m£ PROGRAM OETERniNES TH£ MODE (REAL. INTEGER OR ALPHA) OF EACH WORO C 
ON Th£ TAPE. C 

Th£ program then PRINTS EACH RECORD in the CORRECT FORMAT. c 

c 




DIMENSION NN(S13).FH(3> *FrT (10) * ' ** ** 

OATA FM /7h» A lO » 5 a • 6H»Ei5.7 * 7 h*U0»SX ✓ 

DATA FhT (1) »FMT (10) .CPAREn /AHUIX » iH» • iH> / 

REWIND I 
WRITE (6.1) 

1 rOWMATUHI) 

NRtC=0 

2 1CNT = IVAR(UNN. 0*513) 

IF (ICNT.EO.O) STOP 
IFUCNT.LT. 0) GO TO 7 
NRtC=NREC* l 
wRITE(6,3) nREC 

3 FORMAT UH .«RECORO»U) 

DO 6 I-I.1CNT.B 

DO A Jal.8 
J2=I*J-1 

CAuu WHAT(NN(j2).NTYP£) 

FMTU*J)bFM(NTYPEI 
IF(J2.E0.1CNT> GO TO S 
A CONTINUE 
GO TO 6 

5 FMT(2»J)»CPAREN 

6 wRi 7E(6»fHT) (NN(J)»J»1*J2) 

60 TO 2 

7 WRITE (6*8) 

6 FORMAT (IH *30H****P A R X T Y 6 R R 0 R»**»> 

slop 

END 
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APPENDIX D - CONTINUED 


SUBROUTINE WHATIN.NTYPE) 

C This SUBROUTINE DETERMINES the mode of n 
C NTYP£=1 IF N is ALPHA 

C NTyPE=2 IF N IS real 

C NTYPE-3 if n IS INTEGER 

integer Gd 
DATA LD/AmOOOO/ 
call t-SnFT(BB.-24) 

CALL RSHFT (BB» 24) 

IFCIABSI N > .LE. 9999999999) GO TO 140 
M=N 

CALL RSHFT < M »-24) 

CALL RSHFT < M t 24) 

IF C M .NE.BB) GO TO 100 
M=N 

CALL RSHFT ( M *S4) 

IF ( M .LT.l.OR. M .GT-.47) GO TO 100 

M=N 

CALL RSHFT ( M »-6) 

CALL RSHFT ( H »S4) 

IF t M .LT.l.OR. H .GT.47) GO TO 100 
C ALPHA l 

NTYPE*1 
GO TO A 

C REAL 2 

100 nTYP£=2 

. 00 TO 4 

C INTEGER 3 

140 NTYPE-3 

4 RETURN 

ENO 
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Figure 1. - User Tape Interrogation Using Standard Tape Dump 
Program and Writing to an E Format 
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Figure 1. - Concluded 




Figure 2. - User Tape Interrogation Using Tape Dump Program of Appendix D 





Figure 3. * Viking Dynamic Simulator 
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Computed Input Acceleration 
at Spacecraft Adaptor Interface 



Figure 5. - Viking Configuration Used to 
Determine Input Acceleration -Time 
Histories 
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Figure 6. - Translational Accelerations Used As Input 
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Figure 7. - Typical load Versus Time for a Member of the 
Viking Spacecraft Adaptor Truss 
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Figure 8. - Output Acceleration at the Top of the Viking 
Spacecraft 
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